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We present a semiclassical description of the SU{N) Yang-Mills the- 
ory whose partition function at nonzero temperatures is approximated by 
A\. that of an ensemble of N kinds of interacting dyons. The ensemble is 

^ ' mathematically described by an exactly solvable 3d quantum field theory, 

<^ , allowing calculation of correlations functions relevant to confinement. We 

show that all known criteria of confinement are satisfied in this semiclassi- 
cal approximation: (i) the average Polyakov line is zero below some critical 
^ . temperature, and nonzero above it, (ii) static quarks in any nonzero N- 

fS| ' ality representation have linear rising potential energy, (iii) the average 

f^ , spatial Wilson loop falls off exponentially with the area, (iv) A''^ gluons 

0^ ■ are canceled out from the spectrum, (v) the critical temperature is in good 

^^ I agreement with lattice data. 
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1. Philosophy 

Quantum Chromodynamics (QCD) is hardly an exactly solvable quan- 
tum field theory, even in the large N limit. Therefore, one can either do 
exact calculations in a theory that has more symmetries but is not our world, 
or work with QCD but make approximations. The first is useful as a theo- 
retical laboratory, the second is necessary to understand semi-quantitatively 
the key phenomena, to explain experimental data, and to make predictions. 

An approximation is considered to be legitimate if there is a systematic 
way of improving its accuracy. The semiclassical approach belongs to this 
category. One chooses a saddle-point classical field and then has to take into 
account quantum fluctuations about it. Part of the fluctuations are ultra- 
violet and are thus the same as in empty space. Therefore their role is to 
renormalize the bare coupling constant; at this point the famous dimensional 
transmutation occurs, when the ultraviolet cutoff in a proper combination 
with the bare coupling constant forms the QCD scale parameter A, the 
only dimensional scale that henceforth will be in the theory. What is left, 
is a series in the 't Hooft running coupling A = Nag /27r coming from loop 
expansion in the background of classical configurations. 

The argument of the running coupling A is determined by the largest 
scale in the vacuum, max{T, n*), where T is temperature, and n is the 
mean 4d density of the (large) classical field configurations. For example, 
near the deconfinement temperature T ~ Tc ~ A the running coupling is 
approximately A = -3- In (AirT/e'^^ ) ~ 7 [IJ • The numerically large factor 
27r in the argument of the logarithm is not accidental but related to the 
fact that it is actually not the temperature itself but rather the Matsubara 
frequency 2ttT that defines the scale. At zero temperature many QCD 
specialists believe that Os does not grow above the value of 0.5, which gives 
A ~ ^. Therefore, in the whole range of temperatures within the confining 
phase the semiclassical approximation is expected to yield the accuracy of 
15-25%, already in the 1-loop approximation (provided the saddle point 
is chosen correctly!) with a possibility for rapid improvement when higher 
loops are taken into account. We shall see, however, that the actual accuracy 
can be much better than this estimate. It is not a too big price to pay for 
solving the most challenging riddle in 35 years: confinement. 

We shall be considering the pure Yang-Mills theory based on the SU{N) 
gauge group in a broad range of temperatures between and Tc, the decon- 
finement phase transition temperature. Although the formalism we use is 
designed for nonzero T, we shall see that the physical observables we find 
(such as the string tension) have a finite limit when T — > 0. In this limit 
the nonzero temperature can be thought of as an infrared regulator. After 
all, our world's temperature is 2.7 K 7^ 0. 



Confinement, as we understand it today and learn from lattice experi- 
ments with a pure glue theory, has in fact many facets, and all have to be 
explained. Let us enumerate the main: 

• the average Polyakov line in any A^-ality nonzero representation of the 
SU{N) group is zero below Tc and nonzero above it 

• the potential energy of two static colour sources (defined through the 
correlation function of two Polyakov lines) asymptotically rises linearly 
with the separation; the slope called the string tension depends only 
on the A^-ality of the sources 

• the average of the spatial Wilson loop decays exponentially with the 
area spanning the contour; at vanishing temperatures the spatial ( "mag- 
netic") string tension has to coincide with the "electric" one, for all 
representations 

• the mass gap: no massless gluons left in the spectrum. 

These properties have been obtained in Ref. [2] with Victor Petrov, 
which is the base for this presentation. 

2. Yang Mills theory at nonzero temperatures 

The Yang-Mills (YM) partition function can be written as a path or 
functional integral over the spatial components of the connection Ai(t,x) 
satisfying the periodic boundary conditions up to a gauge transformation 
J7(x) over which one has to integrate separately [3j: 
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where Bi = Cijk [djAk — ^[AjAk]) is the magnetic field strength and A^ 

O^^jO + iQ^di^ is the gauge-transformed potential. Ai, Bi are N xN ma- 
trices belonging to the su{N) algebra while 0, is an element of the SU{N) 
group. 

One can rewrite the partition function in a more customary form by 
introducing gauge-transformed integration variables Ai that are strictly pe- 
riodic in time, and trading 0, for the time component of the YM potential 



A4 that is also periodic: 

Z = JDA^it,^) exp i-^j^dt jd^i^ TrF^^F^J , (2) 

A^[t + ^,^ =^^(i,x), 

where F^y is the usual field strength. This form stresses the fact that 
Euclidean 0(4) symmetry is restored as T ^ 0. 

An important variable is the Polyakov loop: in the formulation ^ it is 
the path-ordered exponent 

L(x) =Pexp(i /"^(it^4(i,x) J (=^(x)). (3) 

In the formulation ([1]) it is nothing but the SU{N) matrix Q{x) over which 
there is a final integration in Eq. ([1]). The eigenvalues of L{x) are gauge 
invariant; we parameterize them as 

L = diag(e2-*'^^e2-*'^^...,e2-*'^^), /ii + . . .+ /x^ = 0, (4) 



and assume that the phases of these eigenvalues are ordered: /Ui < ^2 < 
• • • < fJ'N < fJ'N+i = Mi + 1- We shall call the set of A^ phases {fim} the 
"holonomy" for short. Apparently, shifting //'s by integers does not change 
the eigenvalues, hence all quantities have to be periodic in all //'s. 

The holonomy is said to be "trivial" if L belongs to one of the A^ elements 
of the group center Zn. For example, in SU{3) the three trivial holonomies 
are 



1. /ii = ^2 = /^3 = 



^'-~3'^'-3'^'-3 



f.i---,^2---,f^3-- 




Trivial holonomy corresponds to equal ^u's, modulo unity. Out of all possible 
combinations of /x's a distinguished role is played by equidistant /u's: 

rfr' = -i~ + ^, TVL = 0^ (5) 



For example, in SU{3) it is 
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We shall call it "maximally non-trivial" or "confining" holonomy as it cor- 
responds to TrL = which the 1^* confinement requirement. 

Immediately, an interesting question arises: Imagine we take the YM 
partition function, be it in form ([I]) or ([2]), and integrate out all degrees of 
freedom except the eigenvalues {/Xm} of the Polyakov loop L(x) (or ^^(x)), 
which, in addition, we take slowly varying in space. What is the effective 
action for ^u's? What set of /x's is preferred dynamically by the YM system 
of fields? 

In general, it is a difficult calculational problem that can be addressed 
using various approximations but in one case the result is known exactly. It 
is the case of the supersymmetric M = 1 version of the YM theory (SYM) 
where in addition to gluons there are gluinos in the adjoint representation. 
In order not to spoil supersymmetry one takes not the real temperature 
but rather a Ad space compactified in the time direction, R^ x S^. The 
difference is that in the "real temperature" case one uses periodic conditions 
in the Euclidean time direction for boson fields (gluons) and aniiperiodic 
conditions for fermion fields (gluinos) - that spoils supersymmetry; in the 
"compactification" case one implies periodic conditions for both kinds of 
fields, which supports supersymmetry. However, we shall anyway call the 
inverse circumference of the compactified time direction "temperature" for 
short. 

There is no perturbative contribution to the potential energy in question 
as function of /i's (directly related in this case to the holomorphic super- 
potential) because of the supersymmetric cancelation between boson and 
fermion loops, and the only contribution is nonperturbative coming from 
dyons. It can be reliably computed in the limit of high "temperatures" and 
then claimed to be actually independent of temperature owing to the holo- 
morphy typical in supersymmetry. The result [1] is that the potential energy 
of the system has the minimum at precisely the "maximally non-trivial" or 
"confining" holonomy (0). 

In the non-super symmetric pure YM theory, there is a perturbative ef- 
fective action for slowly varying ^'s. It can be understood as gluon loop(s) 
in the background of a slowly varying field ^4(x). The effective action 
can be expanded in the number of gradients of ^'s. The zero-order term, 
the potential energy with no derivatives, has been computed long ago in 
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Since the piece with no derivatives imphes constant /i's, it has to be propor- 
tional to the 3-volunie V, and hence to T^ by dimensions. pP^''* has exactly 
N zero minima when all /i's are equal modulo unity. Hence, P^"^^^ says that 
at high temperatures the system prefers one of the N trivial holonomies cor- 
responding to the Polyakov loop being one of the A^ elements of the center 
Zj\[, see Fig. 1, top. 




Fig. 1. The perturbative (top) vs. nonperturbative (bottom) potential energy as 
function of the Polyakov line for the SU{2) {left) and SU{2>) {right) groups. The 
perturbative potential energy has minima where the Polyakov loop is one of the 
TV elements of the center Zjv and is maximal at the "confining" holonomy. The 
nonperturbative potential energy, on the contrary, has a single and non-degenerate 
minimum at the confining holonomy corresponding to Tr L = 0. 



However, gradient terms in the effective action indicate that there is a 
problem with the trivial-holonomy points, already at the perturbative level. 



Indeed, the two-derivative term is |T] 
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Hiu) = [i;iu) + V(l -1^)+ 27i?]^od 1 > V'(^) = ^ lnr(z.). (8) 

Since V'(^) ~ — 1/z^ at small z^, the gradient term becomes negative near 
"trivial" holonomy, which signals its instability even in perturbation theory. 
We shall show below that a semiclassical configuration - an ensemble 
of dyons with quantum fluctuations about it - generates a nonperturba- 
tive free energy shown in Fig. 1, bottom. It has the opposite behaviour of 
the perturbative potential energy, having the minimum at the equidistant 
(confining) values of the /i's. There is a fight between the perturbative and 
nonperturbative contributions to the free energy [6j. Since the perturbative 
contribution is ~ T^ with respect to the nonperturbative one, it certainly 
wins when temperatures are high enough, and the system is then forced 
into one of the A^ vacua thus breaking spontaneously the Zjsr symmetry. At 
low temperatures the nonperturbative contribution prevails forcing the sys- 
tem into the confining vacuum. This is the mechanism of the confinement- 
deconfinement phase transition. It is of the second order for N = 2 but first 
order for A'^ = 3 and higher, in agreement with lattice findings. 

3. Dyon saddle points 

Dyons or Bogomolny-Prasad-Sommerfield (BPS) monopoles [7] are (anti) 
self-dual solutions of the nonlinear Maxwell equations, D'^^FJ^^ = 0. In 
SU{N) there are exactly N kinds of 'fundamental' dyons with Coulomb 
asymptotics for both electric and magnetic fields (hence the term "dyon"): 

i^^i /I 0\/00 0\/-100 

±E = b'"^^^x 0-10,01 0,0 00 

^'^1 \ o/voo-i/vo 1 

(9) 
Dyon solutions are labeled by the holonomy or the set of ;u's at spatial 
infinity: 

/ /^i \ 

^4(|xHcx)) ^27rr /i2 (10) 

V /X3 y 

(we illustrate it for the case of SU{3)). The explicit expressions for the solu- 
tions in various gauges can be found e.g. in the Appendix of Ref. [8]. Inside 
the cores which are of the size ~ l/(T(;Um+i — fJ-m)), the fields are large. 



nonlinearity is essential. The action density is time-independent everywhere 
and is proportional to the temperature. Isolated dyons are thus 3d objects 
but with finite action S'dyon = (8vr^/5f^)(/Xm+i — fJ-m) independent of temper- 
ature (here fiN+i = Mi -|- 1). The full action of all A'' kinds of well-separated 
dyons together is that of one standard instanton: Sinst = Svr^/^^. 

In the semiclassical approach, one has first of all to find the statistical 
weight with which a given classical configuration enters the partition func- 
tion. It is given by exp(— Action), times the determinant"^'^ from small 
quantum oscillations about the saddle point. For an isolated dyon as a 
saddle-point configuration, this factor diverges linearly in the infrared re- 
gion owing to the slow Coulomb decrease of the dyon field ([9]) . It means that 
isolated dyons are not acceptable as saddle points: they have zero weight, 
despite finite classical action. However, one may look for classical solutions 
that are superpositions of N fundamental dyons, with zero net magnetic 
charge. The small-oscillation determinant must be infrared-finite for such 
classical solutions, if they exist. 

4. Instantons w^ith non-trivial holonomy 

Remarkably, the needed classical solution has been found a decade ago 
by Kraan and van Baal ^9j and independently and simultaneously by Lee and 
Lu [lUj , see also |llj . We shall call them for short the "KvBLL instantons" ; 
an alternative name is "calorons with nontrivial holonomy". The solution 
was first found for the SU{2) group but soon generalized to the arbitrary 
SU{N) [12j. A nice overview of the solutions has been presented by Pierre 
van Baal at the 2003 School in Zakopane [13]. We shall mention only the 
essentials here. 

The general solution depends on Euclidean time t and space x and is 
parameterized by 3N positions of N kinds of 'constituent' dyons in space 
xi, . . . , xtv and their U{1) phases ipi, . . . , Va^ • All in all, there are AN col- 
lective coordinates characterizing the solution (called the moduli space), of 
which the action S'inst = Svr^/g'^ is in fact independent, as it should be for a 
general solution with a unity topological charge. The solution also depends 
explicitly on temperature T and on the holonomy /xi, . . . , fj,]\f: 

^KvBLL ^ ^a^^^^. ^^^ . . . ^ xjy, V'l, • • ■,^N; T,/Xi, . . . ,/X^). (11) 

The solution is a relatively simple expression given by elementary functions. 
If the holonomy is trivial (all /i's are equal modulo unity) the expression 
takes the form of the strictly periodic 0(3) symmetric caloron [T^ reducing 
further to the standard 0(4) symmetric BPST instanton ^ISj in the T — > 
limit. At small temperatures but arbitrary holonomy, the KvBLL instanton 
also has only a small 0{T) difference with the standard instanton. 
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One can plot the action density of the KvBLL instanton in various cor- 
ners of the parameter (moduh) space, see Fig. 2. 





-5>v^ 



Fig. 2. Action density inside the SU{3) KvBLL instanton as function of time 
and one space coordinate, for large (left), intermediate (middle) and small (right) 
separations between the three constituent dyons. 

When all dyons are far apart one observes A^ static (i.e. time-independent) 
objects, the isolated dyons. As they merge, the configuration is not static 
anymore, it becomes a process in time. In the limiting case of a complete 
merger, the configuration becomes a 4d lump resembling the standard in- 
stanton. The full (integrated) action is exactly the same 5inst = Svr^/^r^ for 
any choice of the dyon separations. It means that classically dyons do not 
interact. However, they do experience a peculiar interaction at the quantum 
level to which we proceed. 

5. Quantum weight of a neutral cluster of N dyons 

Remarkably, the small-oscillation determinant about the KvBLL instan- 
ton can be computed exactly; this has been first done for the SU{2) group in 
Ref. |16j and later generalized to SU{N) in Ref. [T7j. The quantum weight 
of the KvBLL instanton can be schematically written as an integral over 
3A^ coordinates of dyons (the weight does not depend on the U(l) angles 
■0m 1 hence they can be integrated out): 

(12) 
where g is the full 4:N x AN metric tensor of the moduli space, defined as 
the zero modes overlap matrix, and Det(— A) is the functional determinant 
over nonzero modes, normalized to the free one and regularized by the 
background Pauli-Villars method; fi is the Pauli-Villars ultra-violet cutoff 
and g'^ is the bare coupling constant defined at that cutoff. The Jacobian 
det g turns out to be a square of the determinant of an A^ x A^ matrix G^^' 
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such that ^/detg = det G^^' where 

AT 






J |Xjyj Xjn+l| -t |X,jj X,n_i| 



l^m+1 P"rm / , ^m — ^^ 



in=l 



is a matrix whose entries are Coulomb interactions between dyons that are 
nearest neighbours in kind. The Coulomb interactions in the zero mode 
overlap matrix arise naturally from the Coulomb asymptotics of the dyon 
field Q, so it is quite simple to check that Eq. (fT3]) is correct at large sepa- 
rations. A nontrivial fact is that Eq. (J13p is actually exact for all separations 
between dyons, including the case when they strongly overlap like in Fig. 2, 
right. This has been first conjectured by Lee, Weinberg and Yi [18j and 
then proved to be indeed exact at all separations by a direct calculation 
by Kraan [19] and later checked in Ref. [20]. In the last paper it has been 
also shown that in the limit of trivial holonomy (/im = 0) or vanishing tem- 
perature the measure given by Eqs. (|12|13p reduces to the standard instan- 
ton measure written in terms of the conventional "center-size-orientation" , 
which is a rather nontrivial but gratifying statement. 

The functional determinant over nonzero modes Det^ (~^) together 
with the classical action and the Pauli-Villars cutoff combine into the renor- 
malized scale parameter Apv, times a function of dyon separations, Apv and 
T |161ll7j. It is a complicated function which, for the time being, we approx- 
imate by its most essential part: a constant equal to exp (— PP*''^*), where 
ppert jg ^j^g perturbative gluon loop ([7|) in the background of a constant 
field A4 (jlOp . This part is necessarily present in Det~ (—A) as most of 
the 3d space outside the instanton's core is just a constant A4 background, 
and indeed the calculation [161 E] exhibits this piece which is the only one 
proportional to the 3-volume. 

Therefore, we write the weight of the KvBLL instanton i.e. a neutral 
cluster of A^ different-kind dyons as 

Wi ^ Jd^i...d^N detGi /^ exp (-pP-t(/,i, . . . ,;x^)) (14) 

where the fugacity / is 

f = ^^ = 0{N'). (15) 

The bare coupling constant g^ in the pre-exponent is renormalized and 
starts to "run" only at the 2-loop level not considered here. Eventually, its 
argument will be the largest scale in the vacuum, be it the temperature or 
the equilibrium density of dyons. 
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6. Quantum weight of many dyons 

In the vacuum problem, one needs to use not one but 0{V) number of 
KvBLL instantons as the saddle point. Solutions with the topological charge 
greater than 2 will be hardly ever known explicitly as their construction 
runs into the problem of resolving the nonlinear Atiyah-Drinfeld-Hitchin- 
Manin-Nahm constraints. At present 2-instanton solutions characterized 
by a nontrivial holonomy have been found [21] but it is insufficient. Nev- 
ertheless, the moduli space measure of an arbitrary number K of KvBLL 
instantons can be constructed despite the lack of explicit solutions, at least 
in the approximation which seems to be relevant for the large- volume ther- 
modynamics, if not exactly. 

When one takes a configuration of K instantons each made of A^ different- 
kind dyons one encounters also same-kind dyons for which the metric ()13p 
is inapplicable. However, the case of identical dyons has been considered 
separately by Gibbons and Manton [22]. The integration measure for K 
identical dyons following from that work is, in our notations, 

^'''"* = -^y'dxi...dx;^detG|-),, 

(-.ident ^ I 4vrz^m - J2k^i r|x,-Xfc|' ^ = -?' /-j^gN 



T|xi-Xj|' 



iy^j 



where the identity factorial is inserted to avoid counting same configurations 
more than once. 

As in the case of different-kind dyons, this result for the metric can be 
easily obtained at large separations from considering the asymptotics of the 
zero modes' overlap. However, in contrast to the different-kind dyons, it 
is not possible to prove that this expression is correct at all separations. 
Moreover, such an extension of Eq. (J16p is probably wrong. The metric for 
two same-kind dyons has been found exactly at all separations by Atiyah 
and Hitchin |23j : it is more complicated than what follows from Eq. (J16p at 
K = 2 but differs from it by terms that are exponentially small at large sep- 
arations [23]. We shall neglect the difference and use the Gibbons-Manton 
metric at face value. The point is, Eq. (|16p imposes very strong repulsion 
between same- kind dyons (as does the exact Atiyah-Hitchin metric), hence 
the range of the moduli space where the two metrics differ is, statistically, 
not frequently visited by dyons. We do not have a proof that all ther- 
modynamic quantities will be computed correctly with this simplification: 
proving or disproving it is an interesting and important problem. To remain 
on the safe side, one has to admit today that the metric (J16p is applicable 
if the dyon ensemble is sufficiently dilute, that is at high temperatures. 
Nevertheless, physical observables we compute have a smooth limit even at 
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T — > 0. Therefore, it may well prove to be a correct computation at any 
temperatures, but this remains to be seen. 

It is possible to combine the metric tensors for different-kind ()13|) and 
same-kind (|16p dyons into one metric appropriate for the moduli space of 
Ki dyons of kind 1, K2 dyons of kind 2,..., Kn dyons of kind A^. It is a 
matrix whose dimension is the total number of dyons, that is a {Ki -|- . . . -|- 
Kn) X {Ki -I- ... -I- Kn) matrix: 

^mi,nj — OmnUij \^'^^m \^'l 



u ^ \^mi ^m—l,k\ i^ ^ \^mi ^r?i+l,fc| 7^ /■ -^ \^mi ^mk 



5m,n-l Sm,n+1 ^^ 5„ 



-L l^mi ^m+l,j\ ^ \^mi ^m— lj| ^ \^mi ^ 



mi -^m] I 



«7^i 



where Xmi is the coordinate of the i dyon of kind m. Since the statisti- 
cal weight of a configuration of dyons is large when det G is large and small 
when it is small, det G imposes an attraction between dyons that are nearest 
neighbours in kind, and a repulsion between same-kind dyons. The coef- 
ficients -1,2,-1 in front of the Coulomb interactions are actually the scalar 
products of the Cartan generators that determine the asymptotics of the 
dyons' field, see Eq. ([9]). 

The matrix G has the following nice properties: 

• symmetry: Gmi,nj = Gnj,mi 

• overall "neutrality" : the sum of Coulomb interactions in non-diagonal 
entries cancel those on the diagonal: J2nj Gmi,nj = 47ri/m 

• identity loss: dyons of the same kind are indistinguishable, meaning 
mathematically that det G is symmetric under permutation of any pair 
of dyons {i^^ j) of the same kind m. Dyons do not 'know' to which 
instanton they belong to 



• 



factorization: in the geometry when dyons fall into K well separated 
neutral clusters of A'^ dyons of different kinds in each, det G factorizes 
into a product of exact integration measures for K KvBLL instantons, 
detG = (detG(^))'^ where G^^^ is given by eq. ^ 



last but not least, the metric corresponding to G is hyper-Kahler, as 
it should be for the moduli space of a self-dual classical field [23] . In 
fact, it is a severe restriction on the metric. 
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7. Ensemble of dyons 

In the semiclassical approximation we thus replace the YM partition 
function ([2]) by the partition function of an interacting ensemble of an ar- 
bitrary number of dyons of N kinds: 

Z= Y. ^n-^n n/('^--/)detG'(x), (18) 

Ki...Km ^^l-"-'^^'- m=l i=l"' 

where x^j is the coordinate of the i^^ dyon of kind m, the matrix G is given 
by Eq. (fT7|) and the fugacity / is given by Eq. (fT5|) . The overall exponent 
of the perturbative potential energy as function of the holonomy {//} is 
understood, as in Eq. ([HI). 

The ensemble defined by a determinant of a matrix whose dimension is 
the number of particles, is not a usual one. More customary, the interaction 
is given by the Boltzmann factor exp (— [/int(xi, . . .)). Of course, one can 
always present the determinant in that way using the identity det G = 
exp(TrlogG) = exp(— C/int) but the interactions will then include three-, 
four-, five-... body forces. At the same time, it is precisely the determinant 
form of the interaction that makes the statistical physics of dyons an exactly 
solvable problem. 

8. Dyons' free energy: confining liolonomy preferred 

The partition function (fT8|) can be computed directly and exactly, just 
by writing the determinant of G by definition as a sum of permutations 
of products of the matrix entries. The result is astonishingly simple: all 
Coulomb interactions cancel exactly after integration over dyons' positions, 
provided the overall neutrality condition is satisfied, viz. Ki = K2 = . . . = 
Kn = K; otherwise the partition function is divergent. Therefore, the 
recipe for computing the partition function is just to impose the neutrality 
condition and then to throw out all Coulomb interactions! We have thus to 
take the product of (47r^'m)'s from the diagonal of G: 



(4vrW^ - 



The quantity A'KJVis dimensionless and large for large volumes V . The sum 
can be therefore computed from the saddle point in K using the Stirling 
asymptotics for large factorials, and we obtain 

Z = ey.-p{A'KfVN{viV2...UN)'^^ , ui + U2 + . . . + i^N = I- (19) 
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By definition, F = —TlogZ is tlie nonperturbative dyon-induced free en- 
ergy as function of the holonomy; for iV = 2, 3 it is plotted in Fig. 1, bottom. 
Evidently, it has the minimum at 

1 , , 

Ui = U2 = . . . = l^N = JT (20) 

corresponding to equidistant, that is confining values of //'s ([5])! At the 
minimum, the free energy is 

F„,in = -Tlog^^in = -ATTfVT = -^ A^ V = ^ ^ V, (21) 

and there are no Coulomb corrections to this result. In the last equation 
we have introduced the A^-independent 't Hooft coupling A = asN/lir. 

We note that the free energy is 0{N'^) as expected on general A^-counting 
grounds and that it is temperature- independent. It corresponds to \ogZ 
being proportional to the 4- volume V^^' = V/T^ demonstrating the expected 
extensive behaviour at low temperature. 



9. Statistical physics of dyons as a Quantum Field Theory 

Although the Coulomb interactions of dyons cancel exactly in the free 
energy, the dyon ensemble defined by Eq. (jlSp is not a free gas but a highly 
correlated system. To facilitate computing observables through correlation 
functions, we rewrite Eq. (jlSp as an equivalent quantum field theory. As 
a byproduct, we shall also check that the result for the free energy (|19p is 
correct. 

To proceed to the quantum field theory description we use two mathe- 
matical tricks. 

1. "Fermionization" (Berezin [25J). It is helpful to exponentiate the 
Coulomb interactions rather than keeping them in det G. To that end one 
presents the determinant of a matrix as an integral over a finite number of 
anticommuting Grassmann variables: 

det (Gab) = / H '^^l ^^^ ^^P ("^A ^ab i^B 

Now we have the two-body Coulomb interactions in the exponent and 
it is possible to use the second trick presenting Coulomb interactions with 
the help of a functional integral over an auxiliary boson field. 
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2. "Bosonization" (Polyakov [26j). One can write 



exp J2 



^m^r, 



^m -^71 



Dcj) exp 



/^Ki^^^^ 



+ 



expf /p— pj, p = ^Qm<5(x-Xm). 



After applying the first trick the "charges" Qm become Grassmann vari- 
ables but after applying the second one, it becomes easy to integrate them 
out since the square of a Grassmann variable is zero. In fact one needs 2N 
boson fields Vm,Wm to reproduce diagonal elements of G and 2N anticom- 
muting ("ghost") fields xL^Xm to present the non-diagonal elements. The 
chain of identities is accomplished in Ref. [2] and the result for the partition 
function ([18]) is, identically, a path integral defining a quantum field theory 
in 3 dimensions: 

Z = / Dx^ Dx Dv Dw exp d^xl— (dixlndiXm + diVmdiWm) (22) 



+ / 



(-4^P^ + «^)^— + xL g Q 



Xn 



N 



j: _ V^ gUIm-U'm + l 



(23) 



m,=l 



The fields Vm have the meaning of the asymptotic Abelian electric potentials 
of dyons, 



(^4)mn 

^m4(x)/r 



2'KHr. 



2^m(, 



E. 



VA 



mAi 



(24) 



while Wm have the meaning of the dual (or magnetic) Abelian potentials. 
Note that the kinetic energy for the Vm , Wm fields has only the mixing term 
diVmdiWm which is nothing but the Abelian duality transformation E • B. 
The function J-{w) (j23p where one assumes a cyclic summation over m, is 
known as the periodic (or affine) Toda lattice. 

Although the Lagrangian in Eq. (|22p describes a highly nonlinear inter- 
acting quantum field theory, it is in fact exactly solvable! To prove it, one 
observes that the fields Vm enter the Lagrangian only linearly, therefore one 
can integrate them out. It leads to a functional 5-function: 



DVr. 



5 --—d Wm + fj; 

47r OWrry 



(25) 
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This (5-function restricts possible fields Wm over which one still has to inte- 
grate in eq. (j22|) . Let Wm be a solution to the argument of the 5-function. 
Integrating over small fluctuations about w gives the Jacobian 



Jac = det"^ ( -^d'^5mn + f- 




dWrr, dWn 



(26) 



Remarkably, exactly the same functional determinant (but in the numer- 
ator) arises from integrating over the ghost fields, for any background w. 
Therefore, all quantum corrections cancel exactly between the boson and 
ghost fields (a characteristic feature of supersymmetry) , and the ensemble 
of dyons is basically governed by a classical field theory. 

To find the ground state we examine the fields' potential energy being 
—A7rfnmdT/dwm which we prefer to write restoring z/^ = Hm+i — fJ-m and 
J^ as 

r = -AirfV J2 ^m e""^-""^^ (27) 

m 

(the volume factor arises for constant fields Wm)- One has first to find the 
stationary point in Wm for a given set of fm's- It leads to the equations 

8P .„ 



etc. (28) 



dwn. 
whose solution is 

Putting it back into eq. (j27|) we obtain 

V = -ATTfVN{viV2...VN)^^ VI+V2 + ... + VN = l, (29) 

which is exactly what one gets from a direct calculation of the partition 
function, outlined in the previous section, see Eq. ()19p . The minimum is 



achieved at the equidistant, confining value of the holonomy, see Eqs. (l20|5]) . 
Using field-theoretic methods, we have also proven that the result is exact, 
as all potential quantum corrections cancel. It is in line with the exact 
cancelation of the Coulomb interactions in the determinant. 

Given this cancelation, the key finding - that the dyon-induced free en- 
ergy has the minimum at the confining value of holonomy - is trivial. If 
all Coulomb interactions cancel after integration over dyons' positions, the 
weight of a many-dyon configuration is the same as if they were infinitely 
dilute (although they are not). Then the weight, what concerns the holon- 
omy, is proportional to the product of diagonal matrix elements of G in 
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the dilute limit, that is to the normalization integrals for dyon zero modes. 
These are nothing but the field strengths F^u of individual dyons, hence the 
normalization is proportional to the product of the dyon actions ~ Um where 

l^m = /Um+l - ^m and UN = IJ-l + 1 - fJ'N SUch that Z^i + Z/2 + . . . + Z^AT = 1. 

The sum of all N kinds of dyons' actions is fixed and equal to the instanton 
action, however, it is the product of actions that defines the weight. The 
product is maximal when all actions are equal, hence the equidistant or 
confining /i's are statistically preferred. Thus, the average Polyakov line is 
zero, < Tr L >= 0. 



10. Heavy quark potential 

The field-theoretic representation of the dyon ensemble enables one to 
compute various Yang-Mills correlation functions in the semiclassical ap- 
proximation. The key observables relevant to confinement are the correla- 
tion function of two Polyakov lines (defining the heavy quark potential), and 
the average of large Wilson loops. A detailed calculation of these quantities 
is performed in Ref. [2j; here we only present the results and discuss the 
meaning. 



10.1. N-ality and k-strings 

From the viewpoint of confinement, all irreducible representations of the 
SU{N) group fall into A'" classes: those that appear in the direct product 
of any number of adjoint representations, and those that appear in the 
direct product of any number of adjoint representations with the irreducible 
representation being the rank-fc antisymmetric tensor, k = 1, . . . , N—1. "N- 
ality" is said to be zero in the first case and equal to k in the second. A^- 
ality-zero representations transform trivially under the center of the group 
Zj\f] the rest acquire a phase 2'Kk/N . 

One expects that there is no asymptotic linear potential between static 
color sources in the adjoint representation as such sources are screened by 
gluons. If a representation is found in a direct product of some number 
of adjoint representations and a rank-A; antisymmetric representation, the 
adjoint ones "cancel out" as they can be all screened by an appropriate 
number of gluons. Therefore, from the confinement viewpoint all A^-ality 
= k representations are equivalent and there are only N — 1 string tensions 
o"fc,iV being the coefficients in the asymptotic linear potential for sources 
in the antisymmetric rank-/c representation. They are called "/c-strings" . 
The representation dimension is dk,N = k\(N-k)\ ^^^ ^^^ eigenvalue of the 
quadratic Casimir operator is Ck,N = ^^ k{N — k). 
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The value k = l corresponds to the fundamental representation whereas 
k = N — 1 corresponds to the representation conjugate to the fundamental 
[quarks and anti-quarks]. In general, the rank-(A^ — A;) antisymmetric rep- 
resentation is conjugate to the rank-A; one; it has the same dimension and 
the same string tension, a^^N = '^N-k,N- Therefore, for odd N all string 
tensions appear in equal pairs; for even N, apart from pairs, there is one 
privileged representation with k = ^ which has no pair and is real. The 
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total number of different string tensions is thus 

The behaviour of ak,N as function of k and N is an important is- 
sue as it discriminates between various confinement mechanisms. On gen- 
eral A^-counting grounds one can only infer that at large A^ and k <^ N, 
(yk,N/cri,N = ik/N){l + Oil/N"^)). Important, there should be no 0{N-'^) 
correction [27]. A popular version called "Casimir scaling", according to 
which the string tension is proportional to the Casimir operator for a given 
representation (it stems from an idea that confinement is somehow related 
to the modification of a one-gluon exchange at large distances), does not 
satisfy this restriction. 

10.2. Correlation function of Polyakov lines 

To find the potential energy V^^jv of static "quark" and "antiquark" 
transforming according to the antisymmetric rank-A; representation, one has 
to consider the correlation of Polyakov lines in the appropriate representa- 
tion: 

TrLfc,Af(zi) TrL^^^(z2)^ = const, exp i ^'^ i, ^ j . (30) 

Far away from dyons' cores the field is Abelian and in the field-theoretic 
language of Eq. ()22p is given by Eq. ()24p . Therefore, the Polyakov line in 
the fundamental representation is 

TlL{z)='^Zm, Zm = expi2TTifJ,m- -Vmi^)) ■ (31) 

m=l ^ ^ / 

In the general antisymmetric rank-A; representation 

N 

Tr-^fc,Af(z) = 2^ ZmiZm2---Zmk (32) 

mi<m2<...<Tni^ 

where cyclic summation from 1 to A^ is assumed. 

The average (j30p can be computed from the quantum field theory (j22p . 
Inserting the two Polyakov lines (j32p into Eq. (|22p we observe that the 
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Abelian electric potential Vm enters linearly in the exponent as before. 
Therefore, it can be integrated out, leading to a J-function for the dual 
field Wm, which is now shifted by the source (cf. Eq. (j25p ): 

f f T dJ- i 

DVjn > 11'^ ( ~^^^'^™ + -^^ ■^^{x-2.l){dmmi+ ■■■ + Smmk) 

+ -S{x-Z2)idmni+ ■■■ + Smnk)] ■ (33) 

One has to find the dual field Wm{^) nullifying the argument of this 5- 
function, plug it into the action 

exp (|dx ^Hw)) , (34) 

and sum over all sets {rrii < m2 < ... < w-A;}, {ni < n2 < ... < n/.} with 
the weight exp (27ri(mi + . . . + m^ — ni — ... — ni.)/N). The Jacobian from 
resolving the (5-function again cancels exactly with the determinant arising 
from ghosts. Therefore, the calculation of the correlator (j30p . sketched 
above, is exact. 

At large separations between the sources |zi — Z2I, the fields Wm resolving 
the (^-function are small and one can expand the Toda chain: 

1 dJ^ 

Hw) = Y. 6"''"-"''"+^ -N+-Wrn Mmn Wn, T, ~ A^mn Wn, (35) 

m 2 dWm 

where 



M 



/ 2 -1 ... -1 \ 

-1 2 -1 ... 

0-12-1 ... 

V -1 ... -1 2 / 



(36) 



As apparent from Eq. (|35p. the eigenvalues of M determine the spectrum 
of the dual fields Wm- There is one zero eigenvalue which decouples from 
everywhere, and A^ — 1 nonzero eigenvalues 



TTk^^ 



M^'^ = [2sm—j , k = l,...,N-l. (37) 

Certain orthogonality relation imposes the selection rule: the asymptotics of 
the correlation function of two Polyakov lines in the antisymmetric rank-fc 
representation is determined by precisely the k^^ eigenvalue. We obtain [2] 

TrLfc,^(zi) TV 4,^(z2))'''=°° const, exp (-|zi-Z2|Mv/Aiw) (38) 
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where M is the 'dual photon' mass, 

Comparing it with the definition of the heavy quark potential (j30p we find 
that there is an asymptotically linear potential between static "quarks" in 
any A'^-ality nonzero representation, with the /c-string tension 

ak,N = MT^/m^) = 2MT sin ^ = ^ ^ sin ^. (40) 

This is the so-called 'sine regime': it has been found before in certain su- 
persymmetric theories [28]. Lattice simulations [29] support this regime, 
whereas another lattice study [30] gives somewhat smaller values but within 
two standard deviations from the values following from eq. (j40p . For a gen- 
eral discussion of the sine regime for /c-strings, which is favored from many 
viewpoints, see [27] , 

We see that at large iV and A; <C N, ak,N/o-i,N = {k/N){l + 0{1/N'^)), 
as it should be on general grounds, and that all A;-string tensions have a 
finite limit at zero temperature. 

11. Area law for large Wilson loops 

When dealing with the ensemble of dyons, it is convenient to use a gauge 
where ^4 is diagonal (i.e. Abelian). This necessarily implies Dirac string 
singularities sticking from dyons, which are however gauge artifacts as they 
do not carry any energy. Moreover, the Dirac strings' directions are also 
subject to the freedom of the gauge choice. For example, one can choose 
the gauge in which N dyons belonging to a neutral cluster are connected 
by Dirac strings. This choice is, however, not convenient for the ensemble 
as dyons have to loose their "memory" to what particular instanton they 
belong to. The natural gauge is where all Dirac strings of all dyons are 
directed to infinity along some axis, e.g. along the z axis. The dyons' field 
in this gauge is given explicitly in Ref. [8] (for the SU{2) group). 

In this gauge, the magnetic field of dyons beyond their cores is Abelian 
and is a superposition of the Abelian fields of individual dyons. For large 
Wilson loops we are interested in, it is this superposition field of a large 
number of dyons that contributes most as they have a slowly decreasing 
l/|x— Xj| asymptotics, hence the use of the field outside the cores is justified. 
Owing to self-duality, 

T 

[Bi{^)]mn = [diM{^)]^n = "^ ^mn diVm{^), (41) 
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cf. eq. (I24p . Since Ai is Abelian beyond the cores, one can use the Stokes 
theorem for the spatial Wilson loop: 

W = Ti'Pexp i (bAidx' = Tr exp i j Bi cfa' = ^ exp (-i- fcPa' diVmj ■ 

(42) 
Eq. (j42p may look contradictory as we first use Bi = curMj and then Bi = 
diAi^. Actually there is no contradiction as the last equation is true up to 
Dirac string singularities which carry away the magnetic flux. If the Dirac 
string pierces the surface spanning the loop it gives a quantized contribution 
exp(27ri -integer) = 1; one can also use the gauge freedom to direct Dirac 
strings parallel to the loop surface in which case there is no contribution 
from the Dirac strings at all. 

Let us take a flat Wilson loop lying in the {xy) plane at z = 0. Then 
eq. ([12]) is continued as 



W = ^exp[-i— / cfxdzVmSiz) 

= ^expfi^ / cfxvmdjiz)) . (43) 

It means that the average of the Wilson loop in the dyons ensemble is given 
by the partition function (j22p with the source 

Eexp(,f/A.„^.(....A.ea)) 

where 6{x,y G Area) is a step function equal to unity if x, y belong to the 
area inside the loop and equal to zero otherwise. As in the case of the 
Polyakov lines the presence of the Wilson loop shifts the argument of the 6- 
function arising from the integration over the Vm variables, and the average 
Wilson loop in the fundamental representation is given by the equation 

Therefore, one has to solve the non- linear equations on Wm^s with a source 
along the surface of the loop, 

- d^Wm + M2 (e"'--"'-+i - e'"— I-'"'") = -2^i 6^^^ '^^ 0{x, y G Area), 

(45) 
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for all mi, plug it into the action (Anf /N)J^{w), and sum over mi. In 
order to evaluate the average of the Wilson loop in a general antisym- 
metric rank-Zc representation, one has to take the source in eq. (|45p as 
—27ri5'{z) {Smmi+- ■ •+'5mmfc) and sum over rrii < . . . < ruk from 1 to A^, 
see eq. ([32]). Again, the ghost determinant cancels exactly the Jacobian 
from the fluctuations of Wm about the solution, therefore the classical-field 
calculation is exact. 

Contrary to the case of the Polyakov lines, one cannot, generally speak- 
ing, linearize eq. (j45p in Wm but has to solve the non-linear equations as they 
are. The Toda equations (j45p with a S'{z) source in the r.h.s. define "pinned 
soliton" solutions Wm{z) that are Id functions in the direction transverse to 
the surface spanning the Wilson loop but do not depend on the coordinates 
X, y provided they are taken inside the loop. Beyond that surface Wm = 0. 
Along the perimeter of the loop, Wm interpolate between the soliton and 
zero. For large areas, the action (j34p is therefore proportional to the area 
of the surface spanning the loop, which gives the famous area law for the 
average Wilson loop. The coefficient in the area law, the 'magnetic' string 
tension, is found from integrating the action on the solution Wm{z) in the z 
direction. 

The exact solutions of Eq. (j45p for any A^ and any representation k have 

been found in Ref. [2], and the resulting 'magnetic' string tension turns out 

to be 

K^ N Tik 

(^k,N = -. sm -- , (46) 

A vr iV 

which coincides with the 'electric' string tension (j40p found from the corre- 
lators of the Polyakov lines, for all /c-strings! 
Several comments are in order here. 

• The 'electric' and 'magnetic' string tensions should coincide only in 
the limit T ^ where the Euclidean 0(4) symmetry is restored. 
Both calculations have been in fact performed in that limit as we 
have ignored the temperature-dependent perturbative potential ([7]). 
If it is included, the 'electric' and 'magnetic' string tensions split. 



• despite that the theory (j22p is 3-dimensional, with the temperature 
entering just as a parameter in the Lagrangian, it "knows" about the 
restoration of Euclidean 0(4) symmetry at T — > 0. 

• the 'electric' and 'magnetic' string tensions are technically obtained in 
very different ways: the first is related to the mass of the elementary 
excitation of the dual fields Wm, whereas the latter is related to the 
mass of the dual field soliton. 
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12. Cancelation of gluons in the confinement phase 

To prove confinement, it is insufficient to demonstrate the area law for 
large Wilson loops and the zero average for the Polyakov line: it must be 
shown that there are no massless gluons left in the spectrum. We give an 
argument that this indeed happens in the dyon vacuum. 

A manifestation of massless gluons in perturbation theory is the Stefan- 
Boltzmann law for the free energy: 

-|logZsB = ^ = -gTniV^-l). (47) 

It is proportional to the number of gluons A^^ — 1 and has the T^ behaviour 
characteristic of massless particles. In the confinement phase, neither is 
permissible: If only glueballs are left in the spectrum the free energy must 
be 0{N^) and the temperature dependence must be very weak until T tuT^ 
where it abruptly rises owing to the excitation of many glueballs. 

As explained in Section 8, the ensemble of dyons has a nonperturbative 
free energy 

Fa N"^ A^ 

It is 0{N'^) but temperature- independent. We have doubled -Fmin from 
Eq. (j2ip keeping in mind that there are also anti-dyons and assuming that 
their interactions with dyons is not as strong as the interactions between 
dyons and anti-dyons separately, as induced by the determinant measure 
(fT7|) , therefore treating dyons and anti-dyons as two independent "liquids" . 
(By the same logic, the string tension (jlOj) has to be multiplied by \/2 as 
due to anti-dyons.) 

Dyons force the system to have the "maximally nontrivial" holonomy 
([5]). For that holonomy, the perturbative potential energy d?]) is at its max- 
imum equal to 

i^pert.max ^^^4 ( ^2 _ J_\ ^ (49) 

V 45 V N^J 

The full free energy is the sum of the three terms above. 

We see that the leading 0{N'^) term in the Stefan-Boltzmann law is 
canceled by the potential energy precisely at the confining holonomy point 
and nowhere else! In fact it seems to be the only way how 0{N'^) massless 
gluons can be canceled out of the free energy, and the main question shifts 
to why does the system prefer the "maximally nontrivial" holonomy. Dyons 
answer that question. 
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13. Deconfinement phase transition 

As the temperature rises, the perturbative free energy grows as T^ and 
eventually it overcomes the negative nonperturbative free energy (1481) . At 
this point, the trivial holonomy for which both the perturbative and nonper- 
turbative free energy are zero, becomes favourable. Therefore an estimate 
of the critical deconfinement temperature comes from equating the sum of 
Eq. ([IHll and Eq. (^91) to zero, which gives 



N^ A^ 



45 



2^4 Ar4 _ 1 ^2 • 



(50) 



As expected, it is stable in A^. A more robust quantity, both from the 
theoretical and lattice viewpoints, is the ratio Td ^ since in this ratio the 
poorly known parameters A and A cancel out: 



45 



^2Ar2 



x/a Uvr^ (Ar4 



l)sin2^^ 



N- 



1 



vr 



+ 



1 

iV2 



(51) 





5C/(3) 


5C/(4) 


5C/(6) 


5C/(8) 


Tel \fo^ theory 


0.6430 


0.6150 


0.5967 


0.5906 


Tel \fa., lattice 


0.6462(30) 


0.6344(81) 


0.6101(51) 


0.5928(107) 



In the Table, we compare the values from Eq. (j5ip to those measured 
in lattice simulations of the pure SU{N) gauge theories [31]; there is a 
surprisingly good agreement. A detailed study of the thermodynamics of 
the phase transition will be published elsewhere. 



14. Relation to other suggestions to explain confinement 

Several mechanisms of confinement have been suggested in the past. The 
most popular are 

• condensation of monopoles, or the dual Meissner effect [32^ [33] 

• proliferation of center vortices [62\ |31] , see a modern overview 



These two mechanisms and in particular lattice evidence supporting them 
have been reviewed by Jeff Greensite [35] and we are not going to repeat it 
here. What is important, both monopoles and vortices are identified on a 
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lattice by fixing the gauge - choosing the "maximally Abelian" gauge in the 
first case and the "maximally center" gauge in the second. If this gauge- 
fixing procedure is applied to the dyon vacuum of the present paper, the 
maximally Abelian gauge would probably reveal lattice monopoles where 
dyons are placed, and a subsequent application of the maximally center 
gauge would probably reveal center vortices which would be nothing but 
the phantom Dirac strings connecting dyons. Therefore, lattice findings 
that "there is no confinement without Abelian monopoles" and that "there 
is no confinement without center vortices" is presumably in no contradiction 
with the vacuum being formed by dyons. Moreover, recently there have been 
direct observations of dyons on the lattice by the Humboldt Universitat - 
ITEP group, see [36j and further references therein. 

Some time ago we have observed that standard instantons are also capa- 
ble of yielding confinement, provided the instanton size distribution falls off 
as 1/p^ at large sizes p |37j . This regime implies, however, that large-size in- 
stantons inevitably overlap, since in 4d the packing fraction is proportional 
to the fourth moment of the size distribution /j^ which is divergent. There- 
fore, the usual instantons' "center-size-orientation" parameterization being 
all right for dilute systems is inapplicable for the confinement purposes. One 
needs a parametrization of the collective instantons' coordinates that is as 
good for overlapping solutions as it is for dilute ones. 

In an analogous 2d CP^^^ model also possessing instantons such a pa- 
rameterization has long been known: instantons there are parameterized by 
the positions of A^ kinds of "instanton quarks" . The measure of the moduli 
space of multi-instantons is fortunately known exactly [38] and is given by 
a holomorphic function of the instanton quarks' coordinates. The measure 
is invariant under permutation of the instanton quarks (they should not 
'know' what instanton they belong to) and is perfectly valid for overlapping 
instantons, as well as for dilute ones. In the latter case the measure becomes 
the product of instanton "center-size-orientation" measures [39j . 

In the Ad YM theory a similar parameterization of multi-instantons has 
long been sought, starting from the pioneering work of Callan, Dashen and 
Gross who suggested "merons" as instanton constituents [10], but that did 
not work as merons had a divergent action. Zhitnitsky [H], Petrov and my- 
self put much effort in identifying "instanton quarks" for the YM solutions 
but real progress has been achieved in constructing the KvBLL instan- 
tons [9l [To] whose constituents have been found to be the BPS monopoles, 
or dyons. The price is that one is obliged to take nonzero temperatures, 
however if one is interested in the zero-temperature case, T can be con- 
sidered as an infrared regulator which is safe to put to zero at the end, if 
needed. 
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The measure of the multi-instanton space (I17p is now written in terms 
of the coordinates of the constituent dyons. The metric is hyper-Kahler 
(which is the 4d analogue of holomorphy in 2d), the measure is invariant 
under permutation of dyons (they should not 'know' what instanton they 
belong to) and is presumably valid for overlapping instantons, as well as 
for dilute ones. In the latter case the measure becomes the product of the 
instanton "center-size-orientation" measures j2^. Therefore, it seems to be 
the solution of a long-standing problem. 

Two steps in modernizing the semi-classical "instanton liquid" model [6] 
are critical in getting confinement: 

• generalizing instantons in such a way that they can have arbitrary 
holonomy, and allowing nontrivial holonomy, despite that in pertur- 
bation theory it is forbidden 

• writing the quantum weight of instantons with nontrivial holonomy 
through coordinates of constituent dyons, such that it is applicable 
for overlapping instantons. 

What happens, can be summarized as follows: 

• The ensemble of dyons favours dynamically the "maximally nontriv- 
ial" or confining value of the holonomy. This is almost clear, given 
that the weight is proportional to the product of individual actions of 
A'' kinds of dyons 

• Dyons form a sort of Coulomb plasma (but an exactly solvable variant 
of it) with an appearance of the Debye mass both for "electric" and 
"magnetic" (dual) photons. The first give rise to the exponential fall- 
off of the correlation of two Polyakov lines, i.e. to the linear heavy- 
quark potential, the second yield the area law for spatial Wilson loops 

• N^—1 massless gluons cancel out from the free energy, and only massive 
(string?) excitations are left. 

We do not see the quantum- mechanical condensation of monopoles; it is 
hence a new mechanism of confinement. 

15. Why does it work and what should be done next? 

The reason why a semiclassical approximation works well for strong 
interactions (where all dimensionless quantities are, generally speaking, of 
the order of unity) is not altogether clear. A possible justification has been 
outlined in Section 1: After UV renormalization is performed about the 
classical saddle points and the scale parameter A appears as the result of the 
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dimensional transmutation, further quantum corrections to the saddle point 
is a series in the running 't Hooft coupling A whose argument is typically the 
largest scale in the theory, in this case max(T, n^'^) where n is the Ad density 
of the dyons. An estimate shows that the running A is between 1/4 at zero 
temperature and 1/7 or less at critical temperature. Therefore, although 
these numbers are "of the order of unity" , in practical terms they indicate 
that high order loop corrections are not too large. Let us recall that quite 
an accurate computation of anomalous dimensions in critical phenomena 
from the e-expansion by Fisher and Wilson [32] is based on truncating the 
Taylor expansion in e at the first couple of terms, where e = l or sometimes 
20 

Unfortunately, approximations made in Ref. [2] and reproduced above 
are not limited to higher loop corrections. We have (i) ignored dyon inter- 
actions induced by the small oscillation determinant over nonzero modes, 
except the potential energy as function of the holonomy, (ii) ignored the 
interactions of dyons of different duality, treating them as two noninteract- 
ing "liquids", (iii) conjectured a simple form of the dyon measure which 
may be incorrect when two same-kind dyons come close. Although certain 
justification for these approximations can be put forward (see above and 
Ref. [2]) it is desirable not to use them at all, and that may be possible. 

These mathematical problems are of course in the line, as well as fur- 
ther physical problems, probably the most urgent being switching in light 
dynamical quarks into the dyon vacuum, that is moving into the realm of 
the real-world QCD. The main problem there is the spontaneous breaking 
of chiral symmetry. Although we do not think that its mechanism will differ 
dramatically from that found in Ref. [13], as due to the delocalization of 
the near-zero fermion modes, it would be very interesting to see how the 
ensuing effective chiral lagrangian "knows" about the confinement. 
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